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Fragmentation functions are determined for the pion, kaon, and proton by analyzing charged- 
hadron production data in electron-positron annihilation. It is important that uncertainties 
of the determined fragmentation functions are estimated in this analysis. Analysis results 
indicate that gluon and light-quark functions have large uncertainties especially at small . 
We find that next-to-leading-order (NLO) uncertainties are significantly reduced in comparison 
with leading-order (LO) ones in the pion and kaon. The fragmentation functions are very 
different in various analysis groups. However, all the recent functions are roughly within 
the estimated uncertainties, which indicates that they are consistent with each other. We 
provide a code for calculating the fragmentation functions and their uncertainties at a given 
kinematical point of z and by a user. 
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1 Introduction 

A fragmentation function describes a hadronization process from a parent quark, antiquark, 
or gluon to a hadron. Hadron-production processes are often used for investigating important 
physics such as the origin of nucleon spin and properties of quark-hadron matters. Fragmentation 
functions are needed for describing such processes, so that precise functions should be obtained 
for discussing any physics outcome. Nevertheless, it is known that there are large differences in 
the parametrized fragmentation functions, for example, between the ones by Kniehl, Kramer, 
and Potter (KKP)l^and Kretzer^. Recently updated functions by Albino, Kniehl, and Kramer 
(AKK)I3lare also much different from these functions. This fact suggests that the fragmentation 
functions are not determined accurately; therefore, it is important to show reliable regions in 
discussing any hadron-production data. 

Such error analyses have been investigated recently in the studies of unpolarized parton 
distribution functions (PDFs), polarized and nuclear PDFs^^. It is straightforward to apply 
the technique for the fragmentation functions. We determine the fragmentation functions and 
their uncertainties by analyzing the data for charged-hadron production in electron-positron 
annihilation, -\-e~ — > h-\-X. The analyses are done in leading order (LO) and next-to-leading 
order (NLO) of the running coupling constant a^. Because accurate SLD data in 2004 are 



included in our analysis, whereas they are not used in KKP, AKK, and Kretzer's analyses, we 
expect to have improvements. Therefore, important points of our analysis are^ 

• improvement due to addition of accurate SLD data, 

• roles of NLO terms on the determination, namely on the uncertainties, 

• comparison with other analysis results by considering the uncertainties. 

Our analysis method is explained in section 2, results are explained in section 3, and they are 
summarized in section 4. 



2 Analysis method 

The fragmentation function is defined by the ratio of hadron-production cross section to the 
total hadronic cross section: 

fHz,q') = —^^^^^^-^^, (1) 

(Jtot dz 

where is given by the center-of-mass energy squared {Q^ = s), and z is defined by the ratio 
z = Efi/(^/s/2) = 2Ehl ^/Q^ with the hadron energy E^. Since the fragmentation occurs from 
primary quarks, antiquaries, and gluons, the fragmentation function is expressed by the their 
sum: 

F\z, Q2) = a,)D\^zly, Q^). (2) 

Here, Ci{z, Ug) is a coefficient function which is calculated in perturbative QCD, and D^{z, Q^) is 
the fragmentation function of the hadron h from a parton i. The function D^{z, Q"^) is associated 
with a non-perturbative aspect, and it cannot be theoretically calculated in a reliably way. It 
is the purpose of this work to obtain the optimum fragmentation functions for the pion, kaon, 
and proton by analyzing the experimental data for e"*" + e~ h + X. 

In order to determine the functions from the data, we express them in terms of parameters 
at a fixed scale Qq (=1 GeV^): 

D^{z,Ql)=Ntz'^Hl-zf', (3) 

where Nj^, , and (5^ are the parameters to be determined by a analysis of the data. Because 
there is a sum rule due to the energy conservation: Ylh-^t ~ Ylh lo dzzD^{z,Q'^) = 1, it 
is more convenient to choose the parameter M^^ instead of A^^^. They are related by A^^^^ = 
Mj^/B^a'l + 2, /3f + 1), where -B(af + 2,(3^ + 1) is the beta function. In general, a common 
function is assumed for favored functions and different ones are used for disfavored functions. 
The favored indicates a fragmentation from a quark or antiquark which exists in the hadron as a 
constituent in a simple quark model. The disfavored means a fragmentation from other quark or 
antiquark. The details of the formalism are explained in Ref. 5. The optimum parameters are 
determined by minimizing the total given by = J^ji^f"'^"' ~ Fj^'^°)^/((7^"*")^, where Fj^"^**^ 
and F]^^° are experimental and theoretical fragmentation functions, respectively, and 
an experimental error. Uncertainties of the determined fragmentation functions are estimated 
by the Hessian method: 

where Hji^ is the Hessian matrix, is a parameter, ^ indicates the optimum parameter set, and 
the Ax^ value is chosen so that the error becomes the one-a range in the multiparameter space. 
The detailed explanations for the uncertainties are found in Refs. 4 and 5. 



3 Results 



\ Q = Mz 


• SLD 




* ALEPH 




a OPAL 




» DELPHI 




111 



We explain analysis results. First, determined fragmenta- 
tion functions are compared with charged-pion production 
data in Fig. The curve indicates theoretical NLO re- 
sults which are calculated by using determined parameters 
in the analysis, and the uncertainties are shown by the 
shaded band. The comparison suggests that the fit is suc- 
cessful in reproducing the data in four orders of magnitude. 

Determined functions are shown at the initial scales 
(Q^=l GeV^, m^, and ml) and also at an evolved scale 
Q2 = M| in Fig. ^ . The LO and NLO functions and their 
uncertainties are shown. We notice that the uncertainties 
are generally large at small Q^, especially in the LO. The 
gluon and light-quark functions have especially large uncer- 
tainties. However, it is interesting to note that the situation 

is much improved in the NLO because the uncertainties become significantly smaller. The un- 
certainty bands are smaller at large (= M^). Since the fragmentation functions are used at 
small (~ 1 GeV^), for example, in HERMES, RHIC-Spin, and RHIC heavy-ion experiments, 
one should be careful about the reliability of employed functions in one's analysis. 



Figure 1: NLO results are compared 
with charged-pion data. 
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Figure 2: Determined fragmentation functions for tt"*", K'^ , and proton at Q'^ = 1 GeV^, m^, ml, and M|. LO 

and NLO functions are shown with their uncertainties. 



Next, the determined functions are compared with other analysis results for (vr^ + n )/2, 
{K~^ + K~)/2, and {p+p)/2 in Fig. [sP. Our parametrization is denoted HKNS (Hirai, Kumano, 
Nagai, Sudoh). The determined functions in NLO and their uncertainties are shown by the solid 
curves and shaded bands. They are compared with other functions by KKP, AKK, and Kretzer 
at Q^=2, 10, and 100 GeV^. As mentioned earlier, there are much differences between the 
analysis groups. For example, the gluon and s-quark functions have large variations in the pion. 
However, almost all the curves are roughly within the estimated uncertainty bands. It suggests 
that all the analyses should be consistent with each other and that accurate functions cannot 
be determined by the current e^e~ data. After our paper^^, there appeared another analysis by 
de Florian, Sassot, and StratmannEl. Although there are some differences from our functions. 
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Figure 3: Fragmentation functions for (tt^+tt )/2, {K^ + K )/2, and (p + p)/2 are compared with other analysis 

resuhs at Q^=2, 10, and 100 GeV^ 



they are also within the uncertainty bands in Fig. [3l 

The determined fragmentation functions can be calculated by using a code at our web site^ 
by supplying a kinematical condition for z and and a hadron species. It is noteworthy that 
the uncertainties can be also calculated by using the code. 

4 Summary 

The optimum fragmentation functions and their uncertainties have been obtained for the pion, 
kaon, and proton in both LO and NLO of Us by the analyses of charged-hadron production 
data in electron-positron annihilation. It is the first analysis to show the uncertainties in the 
fragmentation functions. The uncertainties were estimated by the Hessian method. We found 
large uncertainties especially at small Q^, so that they need to be taken into account for using 
the functions in the small pT regions of hadron-production measurements in lepton-proton, 
proton-proton, and heavy-ion reactions. We also found that the functions are determined more 
accurately in the NLO than the LO ones particularly in the pion by considering LO and NLO 
uncertainties. There are large differences between previous parametrizations of KKP, AKK, and 
Kretzer, but they are consistent with each other and with our results because they are within 
the uncertainty bands. 
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